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Abstract — Diffusive shock acceleration operating at expanding 
supernova remnant shells is by far the most popular model for 
the origin of galactic cosmic rays. Despite the general consensus 
received by this model, an unambiguous and conclusive proof of 
the supernova remnant hypothesis is still missing. In this context, 
the recent developments in gamma ray astronomy provide us with 
precious insights into the problem of the origin of galactic cosmic 
rays, since production of gamma rays is expected both during 
the acceleration of cosmic rays at supernova remnant shocks and 
during their subsequent propagation in the interstellar medium. 
In particular, the recent detection of a number of supernova 
remnants at TeV energies nicely fits with the model, but it still 
does not constitute a conclusive proof of it, mainly due to the 
difficulty of disentangling the hadronic and leptonic contributions 
to the observed gamma ray emission. In this paper, the most 
relevant cosmic-ray-related results of gamma ray astronomy are 
briefly summarized, and the foreseeable contribution of future 
gamma ray observations to the final solution of the problem of 
cosmic ray origin is discussed. 

I. INTRODUCTION 

Since their discovery in 1912, cosmic ray (CR) nuclei and 
electrons have been studied by means of constantly improving 
direct and indirect detection techniques [1]. Despite exciting 
experimental results and extensive theoretical efforts over the 
past decades, the origin of these particles is still debated. 
The problem is that, unlike photons, CRs are deflected and 
isotropized by the galactic magnetic field and thus their arrival 
direction does not point back to the actual position of their 
accelerators. Revealing the mystery of their origin is of funda- 
mental importance, since CRs can provide unique information 
about the physical conditions of the extreme astrophysical 
objects in which they are likely to be accelerated. Moreover, 
the energy density of CRs, largely dominated by the hadronic 
component, is comparable with the pressure of the galactic 
magnetic field as well as to that of the interstellar medium. 
This implies that CRs have an important role in the dynamical 
balance of our Galaxy [2] and influence interstellar chemistry 
through the heating and ionization of interstellar gas [3]. 

In 1934, Baade and Zwicky first proposed that supernovae 
are the sources of galactic CRs [4]. To support their idea, 
they used a simple argument: the observed CR population can 
be maintained at the present level if a small fraction (a few 
percent) of the galactic supernovae kinetic energy is somehow 
converted into CRs. This argument is strengthened by the 
fact that it is commonly believed that CRs can be efficiently 
accelerated via Fermi mechanism at shock waves that form 



during the expansion of supernova remnants (SNRs) in the 
interstellar medium (see [5]-[8] for reviews). The acceleration 
of CRs in SNRs must be accompanied by a copious gamma 
ray emission due to the decay of neutral pions produced 
in interactions between relativistic protons and protons in 
the interstellar medium. Because the energy transferred to 
accelerated particles is tightly constrained by the observed 
total CR power of the Galaxy, it is possible to obtain an almost 
model independent prediction of the gamma ray luminosity of 
SNRs [9]. 

To date, several SNRs have been detected at TeV energies by 
the major currently operating Cherenkov telescopes (e.g. [10], 
[1 1]-[15] for reviews), whose flux level matches very well the 
above mentioned predictions. Though these results undoubt- 
edly constitute one of the most important advancements in the 
field, they still do not provide us with a definite and direct 
evidence of proton acceleration at SNRs. In fact, competing 
leptonic processes can also explain the observed TeV gamma 
ray emission, provided that the magnetic field does not exceed 
w 10 /iG, and thus accurate modelling is needed in order 
to disentangle the different contributions. Evidence for strong 
w 100 /iG magnetic field, and thus indirect support to the 
hadronic scenario for the gamma ray emission, comes from the 
observation of thin X-ray synchrotron filaments surrounding 
some SNRs [16]-[18] and from the observed rapid variability 
of the synchrotron X-rays from the SNR RXJ1713.7-3946 
[19]. Also neutrinos are produced during the hadronic inter- 
actions responsible for the generation of gamma rays. Their 
detection, though challenging even for the next generation of 
telescopes, would constitute an unambiguous proof for proton 
acceleration in SNRs [9], [20]-[23]. 

Another difficulty of the supernova hypothesis for the origin 
of galactic CRs is related to the very low anisotropy observed 
in the arrival direction of CRs up to energies of lO^^eV 
and above [24]. This is in conflict with the extremely short 
CR confinement time in the Galaxy (which would imply 
anisotropy) required to steepen the CR spectrum from the 
hard injection spectrum (« E^^) predicted by diffusive shock 
acceleration to the steep spectrum (« E"^ '') observed locally 
(for a discussion see e.g. [25]). To solve this problem, not 
only the details of the acceleration mechanism operating at 
SNR shocks needs to be fully understood, but also the way 
in which CRs are released into the interstellar medium during 
the SNR evolution and the way in which they subsequently 



propagate in the Galaxy. 

It is believed that CRs accelerated at SNR shocks are 
gradually released in the interstellar medium as the SNR 
expand and the shock slows down [26] . The bulk of these CRs 
remains diffusively confined within the Galaxy for about ten 
millions years before escaping. During this time CRs undergo 
inelastic interactions in the interstellar gas and produce neutral 
pions which decay into gamma rays (see e.g. [24]). The 
observed galactic diffuse GeV emission which exhibit a good 
spatial correlation with the gaseous disk of the Galaxy is the 
result of these CR interactions [27]. On much smaller scales, a 
correlation between TeV diffuse gamma ray emission and gas 
density has also been observed from regions of the galactic 
disk characterized by the presence of massive molecular cloud 
complexes: the galactic centre ridge [28] and the Cygnus 
region [29]. Such correlations are generally considered as hints 
for a hadronic origin of the gamma ray emission, since the 
presence of massive molecular clouds provide a dense target 
for CR hadronic interactions and thus enhances the expected 
gamma ray emission. A leptonic origin, though not ruled out 
(especially for the Cygnus region [29]), seems disfavored since 
it would require fine tuning to explain the observed correlation. 
For the same reason, the detection of TeV gamma rays from a 
few SNRs spatially associated with dense molecular clouds 
[30]-[33] supports the idea that such TeV emission has a 
hadronic origin and that SNRs might indeed be the sources 
of galactic CRs. 

Even in the absence of a nearby CR accelerator, molecular 
clouds embedded in the "sea" of galactic CRs are expected 
to emit gamma rays. If CRs can freely penetrate the cloud, 
the resulting gamma ray spectrum is expected to mimic the 
underlying CR spectrum. For this reason, molecular clouds 
can be used as "CR barometers" to probe variations of the 
CR spectrum and flux throughout the Galaxy [34], [35]. Such 
variations, if detected, will have to be accounted for by models 
describing the CR injection and propagation in the Galaxy (see 
[36] for a review of propagation models). 

Finally, independently on their nature, the accelerators of 
the galactic CRs are expected to emit gamma rays at some 
level and thus surveys of the galactic plane in gamma rays 
might reveal their presence. Interestingly, more than a half of 
the low galactic latitude sources detected by EGRET at GeV 
energies [37], [38] and many of the TeV sources discovered 
by H.E.S.S. [39], [40] and by MILAGRO [41] along flie 
galactic plane still lack of any clear identification with objects 
at other wavelengths. Revealing the mystery of the nature of 
the unidentified gamma ray sources might shed light on the 
problem of the origin of galactic CRs. 

To conclude, despite the fact that that there is now encourag- 
ing convergence between the model and the observations, we 
are still waiting for the conclusive proof that galactic CRs are 
accelerated at SNRs. The difficulty in distinguishing between 
hadronic and leptonic contribution to the gamma ray emissivity 
of SNRs and the observed isotropy of the diffuse flux of CRs 
up to high energies constitute the two main problems to be 
solved. Forthcoming garmna ray observations in both the GeV 



and TeV domain will constitute a fundamental step towards 
the solution of the problem of the origin of CRs. In particular, 
deep (~ 5 yr) FERMI/GLAST observations of SNRs wiU 
constrain emission models and possibly break the degeneracy 
of leptonic or hadronic origin of the gamma ray emission [42]. 
Moreover, a better determination of the spectrum and spatial 
distribution of the diffuse galactic GeV emission will improve 
our knowledge of the propagation of CRs and their distri- 
bution in the Galaxy. Similarly, TeV instruments of the next 
generation will be required to disentangle the emission form 
leptonic versus hadronic CRs by means of improved spectral 
measurements and morphological studies. An extension of the 
observed energy range up to several hundred TeVs is of crucial 
importance, since at those energies the contribution to the 
observed gamma ray emission from CR electrons is expected 
to be suppressed due to Klein-Nishina effects. Therefore, 
unlike other energy intervals, the interpretation of gamma ray 
observations at these energies is free of confusion and reduces 
to the only possible mechanism: decay of produced in 
hadronic interactions of PeV CRs [22]. 

The paper is structured as follows. In Sec. 2 the available 
gamma ray observations of SNRs and their theoretical interpre- 
tations are reviewed, with particular attention to the question 
on the leptonic or hadronic nature of the emission. In Sec. 3 
the importance of studying molecular clouds in gamma rays 
is highlighted, while Sec. 4 deals with the diffuse gamma ray 
emission observed along the galactic plane and the related 
issue of the spatial distribution of CRs in the Galaxy. In Sec. 5 
the capability of SNRs of accelerating CRs up to PeV energies 
is discussed, and a possible way to test this hypothesis is 
proposed. Conclusions and future perspectives are provided 
in Sec. 6. 

II. GAMMA RAYS FROM SUPERNOVA REMNANTS: 
HADRONIC OR LEPTONIC? 

On average, every « 30 yr a supernova explodes in our 
Galaxy, releasing ~ 10^^ erg in form of mechanical energy 
[43]. If supernovae are the sources of galactic CRs, they have 
to convert « 10% of this energy into relativistic protons in 
order to provide the total power of w 10*^ erg/s needed to 
maintain the galactic CR population at the observed level 
[1], [25]. Such conversion is likely to happen at the SNR 
shocks, where particles are beUeved to be accelerated up to 
ultrarelativistic energies via diffusive shock acceleration [5]- 
[8]. Once the total energy in accelerated protons per SNR 
is determined, it is possible to estimate the expected gamma 
ray luminosity due to hadronic interactions between relativistic 
protons and protons in the intertsellar medium swept up by the 
SNR shock. If we assume a hard CR proton spectrum with 
differential energy distribution Ncr -E^^, as suggested by 
shock acceleration theory, the expected gamma ray emission 
from a single SNR J-y{> Ej) will roughly follow the same 
energy distribution and will depend only on the total energy in 
accelerated protons Wcr, on the density of the ambient gas 



Fig. 1. Left panel: RX J1713.7-3946 as seen by H.E.S.S. (colors) and by ASCA in the lH-3 keV energy band (contours). Figure from Ref. [49]. A shell-like 
morphology is evident at both TeV and X-ray energies. Right panel: spectrum of RX J1713. 7-3945 as measured by H.E.S.S. above 300 GeV. The upper 
solid line is a fit to the data by assuming that gamma rays have a hadronic origin and that the underlying CR proton spectrum is a power law in energy with 
slope 2 and an exponential cutoff at 150 TeV. The lower solid line represent the expected neutrino flux from the same proton population. 



TLgas and on the distance of the SNR d as [9]: 
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These expected fluxes are in general agreement with the fluxes 
of the SNRs observed at TeV energies. 

To date, three young SNRs with apparent shefl-type mor- 
phology have been detected by H.E.S.S. in TeV gamma rays: 
RX J1713.7-3946 (see Fig. 1, left panel) [10], RX J0852.0- 
4622 (Vela Junior) [44] and RCW 86 [45]. These observations 
demonstrate the capability of stereoscopic systems of perform- 
ing detailed morphological studies of extended object at TeV 
energies with an angular resolution of ~ 0.1°. This good 
angular resolution can be exploited to search for similarities in 
the observed morphologies at different wavelengths. Notably, 
for both RX J 17 13.7-3946 and Vela Junior, X-ray observations 
revealed the presence of non-thermal synchrotron emission 
exhibiting a striking morphological similarity with the TeV 
gamma ray image [46], [47]. Such a correlation is naturally 
expected in leptonic models, where both X-rays and gamma 
rays are emitted by the same population of electrons via 
synchrotron and inverse Compton scattering respectively. On 
the other hand, the correlation can be accommodated also 
within hadronic models if most of the emission through both 
7r"-decay gamma rays and synchrotron X-rays comes from 
regions characterized by high magnetic field and gas density 
[1 1], [48]. A scenario in which the X-ray emission is the result 
of synchrotron emission from secondary electrons from decay 



of charged pions would explain the correlation but cannot be 
invoked in this case, since the observed X-ray flux exceeds 
the flux of gamma rays, while one would expect the opposite 
to happen if the electrons were secondary products of proton 
interactions. 

Another remarkable fact is that the gamma ray spectrum of 
RX J1713. 7-3946 has been measured up to a photon energy 
exceeding 30 TeV (see Fig. 1, right panel), providing evidence 
for the existence of primary radiating particles (protons or 
electrons) with spectra extending above such energy [49]. This 
demonstrate the ability of SNRs of accelerating particles up 
to at least w 100 TeV. Such an energy approaches the one of 
the CR knee 4 PeV), where the observed spectrum of CRs 
significantly steepens (see e.g. [24]). Due to the absence of 
any distinct feature in the CR spectrum below the knee, it is 
generally believed that the main contribution to the spectrum 
of galactic CRs comes from a single source population and 
that the representatives of such population must be able to 
accelerate particles up to at least a few PeV. Thus, any 
evidence that SNRs can act as CR PeVatrons will strengthen 
the hypothesis that they indeed are the sources of galactic CRs 
[22]. Future observations in the multi-TeV energy range [50] 
will possibly solve this issue (see Sec. 5 for a discussion). 

Apart from the three objects mentioned above, the list of 
SNRs detected at TeV energies includes also Cas A [51]- 
[53], IC443 [30], [53] and SN 1006 [54]. Good candidates for 
TeV SNRs are HESS J1713-381 (likely to be associated to the 
SNR CTB 37B) [55], HESS 17 14-385 (CTB A) [33] and HESS 
J1745-303 (a possible SNR/molecular cloud association) [32]. 

The question of the hadronic or leptonic nature of the TeV 
emission from SNRs constitutes one of the most discussed 



issues in gamma ray astronomy. If, from the one side, hadronic 
TeV gamma rays are expected if SNRs are the sources of 
galactic CRs, it is also true that most of the SNRs are sources 
of non-thermal X-rays, commonly interpreted as synchrotron 
radiation from multi-TeV electrons (see e.g. [56]-[58]). Such 
electrons can also emit TeV gamma rays via inverse Compton 
scattering off photons of the cosmic microwave background, 
thus providing a competing emission process responsible for 
the TeV radiation (see Fig. 2 for a leptonic model for RX 
J1713.7-3946). If the magnetic field at the SNR is stronger 
than w 10 fxG, then the observed synchrotron X-rays can be 
explained by a relatively meagre number of electrons, which 
would produce unappreciable TeV inverse Compton emission. 
Conversely, if the value of the magnetic field is ~ 10 /zG the 
electrons needed to explain the X-ray emission will also suffice 
to explain the whole observed TeV emission, thus implying 
inefficient acceleration of CR protons. Thus, the value of 
the magnetic field at the shock is a crucial parameter of the 
problem, and its determination would allow us to unveil the 
nature of the gamma ray emission. 

The observation of thin (down to a few arcseconds scale) 
synchrotron X-ray filaments surrounding a number of SNRs 
has been interpreted as an evidence for the presence of a 
strong (« 100 /uG or more) magnetic field at the shock (see 
e.g. [16]-[18], [59]-[62]). According to this interpretation, the 
formation of such filaments is due to the fact that, in such 
a strong magnetic field, the synchrotron cooling time of X- 
ray emitting electrons is very short. Thus electrons radiate X- 
ray synchrotron photons before being significantly advected or 
diffuse away downstream of the shock [63]-L65J. The value 
of the magnetic field is estimated by comparing the observed 
width of the filament with the expected one, which is roughly 
~ VkdTd, where kd and Td are the diffusion coefficient and 
synchrotron loss time of the emitting electrons, both quantities 
depending on the value of the magnetic field [66]. The derived 
value of the magnetic field relies on the (very plausible) 
assumption that diffusion is proceeding at the Bohm rate, 
even though attempts to use X-ray data to derive the diffusion 
coefficient have been made [65], [67]. However, such estimates 
[65], [67] are seriously affected by the uncertainty in the 
determination of supernova parameters such as the shock speed 
and, most of all, the cutoff energy in the synchrotron spectra. 

The presence of a high magnetic field at the shock indirectly 
supports the hadronic origin of the TeV emission for two 
main reasons. First, according to theoretical studies, magnetic 
field amplification due to non-resonant CR driven instability 
is expected to take place at shocks which are accelerating 
efficiently CR protons [68], [69]. The predicted value of the 
amplified field, higher than earlier estimates based on CR 
generation of resonant Alfven modes [70], is compatible with 
the values inferred from X-ray observations and it is also 
the value needed to allow acceleration to proceed up to the 
energy of the CR knee or even above [69]. Second, as said 
above, if the magnetic field is stronger than ~ 10 /zG the 
inverse Compton scattering contribution to the TeV emission 
is negligible, leaving the hadronic channel as the only viable 



mechanism for gamma ray production. 

A high value of the magnetic field is not the only possi- 
ble interpretation for the observed narrow X-ray synchrotron 
filaments. Filaments can also be formed due to damping of 
magnetic turbulence downstream of the shock [71]. In this 
scenario, though magnetic field amplification can still operate 
at the shock , the width of filaments can not be used to estimate 
the value of the magnetic field, instead it simply reflects the 
spatial structure of the magnetic field. An important difference 
between the magnetically Hmited filaments and the energy 
loss limited filaments is that in the former case filamentary 
structures, though broader and of less amphtude, should be 
observed also at lower (e.g. radio) frequencies [71], while in 
the latter case no filaments are expected at radio frequencies 
since low energy electrons are unaffected by synchrotron 
losses [64], [73]. In the case of the SNR SN 1006 the profile 
of the radio emission is quite broad, with only shallow bumps, 
when compared to the narrow peak in the X-ray emission [72]. 
However, an accurate comparison between observation and 
theory is needed in order to draw firm conclusions. Similarly, 
a comparison between the X-ray and radio morphology of 
the Tycho SNR has been attempted, but results stiU remain 
inconclusive [59]. 

Other evidences for the presence of a strong magnetic 
field at SNR shocks comes from the recently discovered fast 
variability of synchrotron X-rays hot spots in the shell of 
RX J17 13. 7-3946 [19] and Cas A [74] . An extremely strong 
magnetic field of « 1 mG has been inferred by comparing 
the decay time of the X-ray flux, of the order of a year, 
with the synchrotron cooling time. Though the estimate of 
the magnetic fields refers only to the sub-parsec scale hot 
spot regions which are seen to vary in X-rays, it suggests 
that significant amplification of the field can possibly happen 
throughout the whole SNR shell. 

Finally, another way to discriminate between hadronic or 
leptonic origin of the TeV emission consists in comparing 
model predictions with the observed X-ray and TeV gamma 
ray spectra. While hadronic models based on predictions 
of non-linear shock acceleration theory generally fit quite 
satisfactorily both X-ray and gamma ray data (for RX J 17 13.7- 
3946 see e.g. [48], [75], [76] and right panel of Fig. 1 
for a fit of H.E.S.S. data based on hadronic interactions), 
leptonic models seem to provide worse fits (see Fig. 2). High 
quality SUZAKU data have been obtained over two decades 
in energy (« 0.5 ^ 50 keV) for die SNR RX J1713.7-3946 
[77]. Such data allow us to constrain the spectrum of the 
synchrotron emitting electrons with unprecedented accuracy 
and over a wide energy range. Remarkably, the resulting 
electron spectrum is very close to that expected fom elec- 
trons accelerated via diffusive shock acceleration in presence 
of synchrotron losses (see Fig. 2) [57], [77]. The electron 
spectrum derived in this way can be used to calculate the TeV 
spectrum due to inverse Compton emission off photons of the 
cosmic microwave background and of the galactic optical and 
infrared background. It is clear from Fig. 2 that such an inverse 
Compton spectrum fails to fit the lowest energy part of the 




Fig. 2. Observed spectram of RX J1713. 7-3946 (data points) with a leptonic model (solid lines). Synchrotron and Inverse Compton (IC) emission are shown 
for an injection spectrum of electrons which is a power law in energy with slope 2. An excellent fit to the X-ray data is obtained by assuming that the shape 
of the cutoff in the electron spectrum is the one expected from diffusive shock acceleration in presence of synchrotron losses [57]. Magnetic field is 14 /tG. 
The inverse Compton flux has been calculated by using the interstellar radiation model from Ref. [78]. The model fails to fit the low energy part of the 
measured gamma ray spectrum. Figure from Ref. [77]. 



measured TeV spectrum, even assuming an optimistic model 
for the interstellar radiation field [78]. The fit can be improved 
only at the expense of enhancing the infrared background of 
a factor of w 20-^ 100 with respect to the adopted value or by 
adding a second electron component. Both these possibilities 
seems quite artificial and thus an hadronic origin of the TeV 
emission seems favored based on spectral information only 
[76], [77]. 

Recently, a number of order of magnitude estimates have 
been presented against hadronic models for the origin of the 
TeV emission from SNRs, and thus indirectly supporting the 
leptonic scenatio [79]-[81]. The main criticisms against the 
hadronic models can be summarized as follows, taking as 
reference case the SNR RX J1713.7-3946: i) in order to fit the 
TeV data within an hadronic model, the gas density cannot 
be too low (see Eq. 1). The high density would then imply 
a strong X-ray thermal Bremsstrahlung emission which is 
not observed [79]. ii) The high value of the magnetic field 
expected in hadronic models implies that a meager number 
of electrons is needed to explain the observed synchrotron 
X-ray emission. This translates into very low values of the 
ratio between CR electrons and protons K^p in the SNR, 
significantly smaller that the value K^p ~ 10^^ observed 
in the Galaxy [79], [80]. Hi) There is no clear tight spatial 
correlation down to small angular scales between the gamma 
ray flux and dense molecular clouds, while this should be 
observed if gamma rays had a hadronic origin [81]. 

The issue of the missing thermal emission (problem i)) 
can be solved by reminding that the shock heating of the 
plasma is relevant only for ions, which carry the inertia of the 



flow. Electrons can be heated downstream of the shock via 
Coulomb collisions with hot ions, but the characteristic time 
scale of the process is much too long to establish electron-ions 
temperature equilibrium [82], [83]. Other processes, possibly 
mediated by plasma waves, might heat electrons and facilitate 
the equipartition (see e.g. [84]), but since the nature of such 
processes is uncertain, it is not possible to draw any firm 
conclusion. It has to be noted that the plasma temperature can 
be further reduced if the shock is converting a large fraction of 
the available kinetic energy into CRs, since this has to happen 
at the expenses of a reduced gas heating [85], [86]. This would 
substantially reduce the thermal X-ray emission, since the peak 
of the emission would be shifted towards the UV energy range 
[87]. For RX J1713. 7-3946, if an electron temperature equal 
to one hundredth of the proton temperature is assumed, no 
thermal X-rays are expected [76]. 

The low value of the electrons to protons ratio K^p (problem 
ii)) is not necessarily a problem for hadronic models. In fact, 
the Kep value of a SNR at a given age does not need to 
agree with the value measured for local CRs. Low energy 
electrons are likely to be accelerated during the late phase 
of the SNR evolution, when the K^p value can be different 
from the present one [77], Also the lack of a good spatial 
correlation between gamma ray emission and the position of 
dense clouds (problem Hi)) can be explained by remembering 
that molecular clouds are mainly neutral. This might reduce 
the level of magnetic turbulence (via damping) needed to have 
effective particle acceleration [76], [88] and would worsen 
the spatial correlation between gamma rays and gas density. 
Independently of the actual nature of the emission, the im- 



proving quality of multiwavelength data and the intrinsic non- 
linearity of the problem of efficient CR acceleration [7], [8] 
make the above mentioned simplistic calculations inadequate 
and demand for a careful and detailed modeling [89]. 

At GeV energies, a number of spatial associations between 
EGRET sources and SNRs have been proposed [90], [91], but 
despite extensive investigations, not a single SNR has been 
unambiguously identified by EGRET. However, a correlation 
was claimed to exist between unidentified EGRET sources and 
SNRs [92]. The detection of the SNR IC443, first proposed 
in Ref. [90] and recently confirmed by the AGILE team 
[93] constitutes the first clear evidence for a GeV SNR. The 
detection of SNRs at GeV energies is of crucial importance 
for constraining emission models. An extended coverage of 
the gamma ray measurements from GeV to TeV energies 
will faciUtate the distinction between hadronic and leptonic 
models, since they predict quite different spectral features (see 
Figs. 1 & 2). In particular, the detection of the slight spectral 
hardening of the spectrum expected in hadronic models [75], 
[89], would allow to test non Unear shock models for CR 
acceleration, which predict concave spectra for the accelerated 
particles [7], [8], [94], [95]. Simulations of FERMI perfor- 
mances indicates that deep observations (~ 5 yr) are needed 
to firmly detect TeV bright remnants such as RX J1713.7-3946 
if the emission is hadronic [42]. 

Finally, the detection of neutrinos would constitute an un- 
ambiguous signature for the presence of CR protons in SNRs. 
However, at least for RX J17 13.7-3946, the neutrino spectrum 
is expected to cut off at an energy of few TeV, if the gamma ray 
emission is entirely due to hadronic interactions (Fig. 1, right 
panel). Unfortunately this is the energy region where neutrino 
telescopes reach their best performances. Thus, a detection of 
SNRs in neutrinos, though possible, seems challenging [21]. 

III. MOLECULAR CLOUDS: PROBES OF THE CR 
SPECTRUM THROUGHOUT THE GALAXY 

The importance of the detection of molecular clouds in 
gamma rays is widely recognized, especially in relation to 
the problem of the origin of CRs. Molecular clouds located 
in the vicinity of CR accelerators, such as for example SNRs, 
provide a dense target for CR hadronic interactions, amplifying 
the resulting gamma ray emission from neutral pion decay 
and making easier the identification of CR sources [96], [97]. 
On the other hand, even in the absence of a nearby CR 
accelerator, molecular clouds are still expected to emit gamma 
rays at a certain level, due to the interaction of background 
CR that penetrate the cloud [34], [35]. According to both these 
scenarios, molecular clouds are passive gamma ray emitters, 
in the sense that they provide a dense target for interactions 
of CRs which are accelerated somewhere else [98]. It has also 
been suggested that particle acceleration can operate inside 
molecular clouds, due to the presence of strong magnetic 
turbulence that might effectively scatter particles [99]. 

The hadronic gamma ray emission from a SNR might be 
enhanced if the supernova shock is overtaking a massive 
molecular cloud [100]. This is because the molecular cloud 



provides a thick target for proton-proton interactions and the 
related gamma ray emissivity is expected to scale hnearly with 
the gas density (see Eq. 1). For this reason, the observed spatial 
associations between a few TeV bright SNRs and massive 
molecular clouds [30]-[33] is considered as a fact supporting 
the hadronic origin of the emission and thus also the idea 
that SNRs might be the sources of galactic CRs. However, the 
connection between high gas density and enhanced gamma 
ray emissivity is not as straightforward as it might appear. In 
fact, the gas constituting molecular clouds is mostly neutral, 
which means that the magnetic turbulence on which particles 
acceleration relies, can be effectively damped [101], [102], 
potentially reducing both the acceleration efficiency and the 
maximum energy of accelerated particles. This might have an 
important effect in suppressing the TeV gamma ray emission 
(see e.g. [88]). Thus, a tight spatial correlation between gas 
density and gamma rays is not necessarily to be expected, 
even if the TeV emission from SNRs has a hadronic origin. 
In the future, high resolution gamma ray observations of 
molecular clouds overtaken by SNR shocks wiU be useful to 
assess the relative importance of these two crucial effects: the 
enhancement of the gamma ray emission due to high density, 
and its reduction due to damping of magnetic turbulence and 
subsequent suppression of the acceleration efficiency. 

The association between SNRs and molecular clouds can 
also be used to determine the distance of the remnant itself, 
which is in most cases quite uncertain. Conversely, the distance 
of a molecular cloud can be inferred from the measured 
doppler shift of CO Unes, coupled with the knowledge of the 
galactic rotation curve [103]. Notably, for the remnant RX 
J17 13.7-3946, the possible correlation of TeV gamma rays 
with CO data obtained with the NANTEN sub-mm telescope 
1 104], allows to estimate the distance of the object of about 
1 kpc. Correspondingly, an age of 1000 3000 years can be 
deduced, which supports the idea that RX J1713. 7-3946 is the 
remnant of the supemova that, according to Chinese historical 
records, exploded in 393 AD [105]. 

Even if the molecular cloud is not overtaken by a SNR shock 
but it is located at a given distance Id from the remnant, it 
can still be illuminated by CRs that escape from the SNR 
and produce there gamma rays [22], [106]. For this scenario, 
it has been shown that, for typical SNR parameters and for 
a distance d ~ 1 kpc, a molecular cloud of mass W^Mq 
emits TeV garmna rays at a detectable level if it is located 
within few hundred parsecs from the SNR [22]. In this case 
the angular displacement between the remnant and the gamma 
ray emission is « 6°{d/l kpc)~^(ici/100 pc). This translates 
into the fact that sometimes the association between a gamma 
ray bright molecular cloud and the SNR which is accelerating 
the CRs responsible for the emission can be not so obvious, 
given that the separation between the two objects can be 
bigger than the instrument field of view. More massive clouds 
(si IO^Mq or above) irradiated by CRs might be detected up 
to the distance of the Galactic centre or beyond. 

Non thermal emission from radio to X-rays is also expected 
from molecular clouds irradiated by CRs coming from a 
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Fig. 3. Broad band spectrum for a MC of mass 10''Mq, radius 20 pc, density ~ 120 cm"^, magnetic field 20 fiG. The MC is at 100 pc from a SNR that 
exploded 2000 yr ago. Distance is 1 kpc. The dotted line shows the emission from Tr^-decay (curve 3), from both background CRs and CRs escaping from 
the remnant. The dot-dashed lines represent the synchrotron (curve 2) and Bremsstrahlung (curve 4) emission from background CR electrons that penetrate 
the MC and the dashed lines the synchrotron (curve 1) and Bremsstrahlung (curve 5) emission from secondary electrons, respectively. Figure from Ref. [107]. 



nearby SNR, due to synchrotron and Bremsstrahlung radiation 
from secondary electrons produced during the same hadronic 
interactions responsible for the production of gamma ray 
photons [107], [108]. In Fig. 3 (from Ref. [107]), as an 
illustrative example, the broad band emission from a molecular 
cloud is plotted from radio waves to gamma rays. The cloud 
has a mass of lO^M©, a magnetic field of 20 /iG and is located 
100 pc away from a SNR which is converting 30% of its total 
mechanical energy 10^^ erg) into CRs. The age of the 
SNR is 2000 yr and the distance to the observer is 1 kpc. 
The contributions to the total radiation from background CRs, 
CRs escaped from the nearby SNR and secondary products 
of CR hadronic interactions in the cloud are plotted, as 
indicated in the figure caption. It is clear that the gamma 
ray emission by far exceeds the emission in the other energy 
bands. Following this rationale, it was proposed that some of 
the unidentified TeV sources detected so far, with no obvious 
or very weak counterparts in other wavelenghts (the so called 
"dark sources") [39]-[41], might be in fact associated with 
massive molecular clouds illuminated by CRs [22]. Such a 
scenario is supported by the fact that many of the unidentified 
TeV sources are spatially extended, as molecular clouds are 
expected to appear. The same interpretation has also been 
suggested for the TeV gamma ray emission detected from a 
molecular cloud associated with the SNR IC 443 [109]. 

In general, the hadronic gamma ray emission from a cloud 
illuminated by a nearby SNR accelerating CRs is expected 
to be the sum of two contributions: the gamma rays form 
interactions of CRs coming from the nearby SNR and the 
gamma rays from interactions of background CRs. The former 
component is expected to exhibit a hard spectrum, reflecting 
the fact that low energy particles remain confined in the 
remnant for long time and thus cannot reach the cloud, while 
the latter will have a steep spectrum, reflecting the underlying 
proton spectrum of the background CRs with slope ~ £^~^^. 



Thus, the gamma ray emission from a molecular cloud located 
in the proximity of a SNR might be characterized, at least 
at some stage of the remnant evolution, by a very peculiar 
spectrum which is steep at low (GeV) energies and hard at high 
(TeV) energies [107]. This can be clearly seen in Fig. 3, where 
the low energy bump in the gamma ray emission, peaking at w 
100 MeV - 1 GeV energies, is the result of the interaction of 
background CRs, while the high energy bump, extending up to 
TeV energies and beyond, is due to interactions of CRs coming 
from the SNR. The possibility of detecting sources with such a 
spectrum is also relevant for the issue of identifying GeV and 
TeV unidentified sources. One of the criteria generally adopted 
to support an association between a GeV and TeV source is, 
beside the positional coincidence, the spectral compatibility 
[110]. Such a criterium would not be applicable to the case 
described here. 

Let us now consider the case of an isolated molecular cloud, 
with no CR sources located in its proximity. In this case the 
molecular cloud is embedded in the "sea" of Galactic CRs, 
assumed to be constant throughout the Galaxy. If CRs can 
freely penetrate into the cloud, then the total gamma ray flux 
above a given energy is solely determined by the mass of 
the cloud and its distance to the Earth and reads [35]: 

Fcii> E^) ~ 1.45 X 10-"i;-i^75 (^^^ cm-^s-i (2) 

where M5 is the cloud mass in units of lO^M© and dkpc is 
the cloud distance in kpc. This means that, since the mass 
and the distance of a cloud can be determined by means 
of observations of CO lines (see e.g. [103], [111]), gamma 
ray observations of molecular clouds can be used to probe 
variations of the CR spectrum in different regions of the 
Galaxy. The condition for the detectabiUty of passive clouds 
embedded in the CR background with EGRET is 5 Mr, /dl^^ > 
10, where 6 represents an enhancement (or suppression) of 



the CR density with respect to the local value. The EGRET 
diffuse gamma ray emission from the direction of prominent 
and nearby molecular clouds has been analyzed, leading to 
a rough agreement with the expactations from models that 
assume a constant CR spectrum throughout the Galaxy (6 « 1) 
[112]. However, since there are only a few nearby clouds in the 
Galaxy satisfying this condition, a more sensitive instruments 
Uke FERMI is needed for these studies. 

At TeV energies, several molecular clouds have been de- 
tected in the galactic centre region |28J or in the neighbour- 
hood of a few SNRs [30J-L33]. In the cases in which the 
mass and the distance of the cloud could be measured, it has 
been possible to evaluate the spectrum and the normalization 
of the CRs at the position of the cloud, by assuming that 
the emission has an hadronic origin. Such an assumption is 
very well justified, especially in the case of the galactic centre 
ridge, because a good spatial correlation is observed between 
TeV emission and gas density, and because the few-degree- 
scale extension of the emission is hard to be explained as 
emission from electrons, since thay suffer severe synchrotron 
losses and thus cannot propagate large distances away from 
the acceleration regions. For the galactic centre region an 
enhancement in the CR density of a factor of 5 = 3 -j- 9 
with respect to the local CR spectrum has been inferred from 
gamma ray observations [28], while for the molecular clouds 
detected in the field of the SNR W28 the inferred overdensities 
of CRs lies in the range 5 = 10 -j- 30 [31]. This overdensities 
refer to CRs with energy above 1 TeV, which are the ones 
producing the observed gamma ray emission above 100 GeV 
In both cases, the inferred spectrum of the CRs is harder than 
the one observed locally. These results demonstrate that there 
are places in the Galaxy, as the galactic centre ridge, where 
the CR spectrum can significantly differ from the local one. 

More generally, a criterium for the detectability of molecular 
clouds with Cherenkov telescopes can be derived by recalling 
that such objects in general have an angular extension signif- 
icantly larger than the typical instrumental angular resolution 
which is of the order of OpsF ~ 0.1°. For a Cherenkov 
telescope Uke H.E.S.S. the sensitivity at 1 TeV for an extended 
sources of angular size ds after 50 hours of exposure is 
roughly « 10""(6's/6'psf) TeV/cm^/s. A cloud with mass 
10^ M5 Mq and typical (and, for simpUcity here, uniform) 
density 100 cm~^ located at a distance dkpc kpc will 
appear as a very extended gamma ray source with angular 
size 0s ~ 1° m}J / dkpc- By combining this information with 
Eq. 2, and recalling that the quantity 5 represents the excess 
of multi-TeV CRs with respect to the local value, it is possible 
to obtain a condition for the detectability that reads: 

if expressed as a function of the point source sensitivity ^^^^ 
of a generic instrument in units of the point source sensitivity 
of H.E.S.S. (fijESS- This is a very rough order of magnitude 
estimate, but it can still be used as a rule of the thumb to 



assess the capabilities of Cherenkov telescopes in detecting 
extended molecular clouds. 

It is clear from Eq. 3 that the detection of TeV gamma rays 
from massive (M5 « 1) molecular clouds embedded in the 
galactic CR background {5 « 1) is beyond the capabilities of 
currently operational Cherenkov telescopes like H.E.S.S.. This 
is in Une with the fact that the molecular clouds detected so 
far in TeV gamma rays are very massive (M5 >> 1) and/or 
located in regions characterised by an excess in the CR density 
{5 » 1). Next generation instruments, such as CTA [113] 
and AGIS [114] , with an expected sensitivity improved by 
about one order of magnitude will detect more distant and 
weaker molecular clouds, and possibly nearby clouds even 
in the less optimistic case when 5 = 1. This will lead to 
an improved knowledge of the distribution of CRs in the 
Galaxy and will possibly provide stringent constraints for CR 
propagation models. 

Of course, what is discussed above is valid only under the 
assumption that CRs freely penetrate the clouds. The issue 
of the penetration or exclusion of CRs from clouds has been 
investigated in several papers [99], 1 1 15], [1 16], in which quite 
different conclusions have been drawn, going from the almost 
free penetration to the exclusion of CRs up to tens of GeV. 
The exclusion of CRs from clouds can suppress the expected 
gamma ray emission and lead to a spectral hardening towards 
the centre of the cloud, due to the fact that low energy CRs 
are excluded more effectively than high energy ones [117]. 

IV. THE DIFFUSE GAMMA RAY EMISSION 
FROM THE GALACTIC PLANE 

The diffuse GeV emission detected by EGRET along the 
galactic plane is believed to be mainly due to the decay 
of neutral pions produced in hadronic interactions between 
CRs and interstellar matter. A strong point in favor of this 
hypothesis is the fact that the diffuse GeV emission correlates 
with the spatial distribution of the interstellar gas which 
constitutes the target for CR interactions [27]. Below ~ 100 
MeV inverse Compton and electron Bremsstrahlung become 
the main contributors to the gamma ray emission. 

Several attempts to model the diffuse galactic gamma ray 
emission have been developed so far, from simple earher mod- 
els where the CR flux was assumed to be constant throughout 
the Galaxy [118] or determined by using dynamical balance 
arguments [119], to numerical studies of CR propagation in 
the Galaxy, tuned in order to match the CR spectrum observed 
locally 1 120]. Remarkably, all these models reproduce quite 
well both the level and the spectrum of the diffuse emission 
with the notable exception of the so called "GeV excess", 
where data exceed model predictions by roughly a factor of 2. 
The GeV excess can be explained by relaxing the assumption 
that the spectrum of CRs observed locally (nuclei, electrons 
or both) is representative for the whole Galaxy. Data can be 
indeed fitted by assuming that the actual CR proton spectrum 
in the Galaxy is harder than the one measured at Earth [121], 
[122]. Another possibility is to assume a hard spectrum of 
CR electrons or that the typical CR electron intensity in the 



Galaxy is significantly higher than that measured locally [120]. 
This would make the inverse Compton contribution to the total 
diffuse gamma ray emission dominate over the hadronic one 
at photon energies above 10 100 GeV. 

Recently, the MILAGRO collaboration reported on the 
detection of multi-TeV diffuse gamma-ray emission from two 
regions in the Galactic plane: the Cygnus region (located 
at galactic longitude 65° < Z < 85°) [29] and the portion 
of the inner Galaxy visible from the location of MILAGRO 
(30° < I < 65°) [123]. When placed alongside the above 
mentioned lower energy measurements by EGRET, such de- 
tections provide an interesting insight into both the spatial 
distributions of CRs in the Galaxy and their spectrum over a 
broad energy range. PuzzUngly, conventional CR propagation 
models that reproduce quite fairly at least the observed level 
of the GeV emission (with the exception of the GeV-excess), 
fail to reproduce the MILAGRO observations by about one 
order of magnitude in the case of the Cygnus region [29] and 
a factor of ^ 5 for the inner Galaxy [123]. However, it has 
to be kept in mind that such models are designed to describe 
the global gamma-ray emission from the Galaxy, and are not 
expected to reproduce correctly local excesses in the gamma 
ray emission due, for example, to the presence of locaUzed 
CR sources. Thus, the discrepancy between models predictions 
and MILAGRO data is representative of the fact that the CR 
spectrum measured locally is probably not representative for 
the whole Galaxy. As done in the case of the GeV excess, a 
harder spectrum of CR protons or electrons can be invoked to 
solve the problem. 

Thus, the origin of the multi-TeV emission detected by 
MILAGRO from the Cygnus region and from the central 
Galactic plane region remains unclear. Either leptonic or 
hadronic processes have been proposed to be the dominant 
contributors to the detected gamma-ray flux, with the possi- 
bihty of a transition from a dominant hadronic contribution 
to a leptonic one in the GeV-TeV energy range [120]. In 
the case of the hadronic scenario, one would expect, as is 
observed, a narrow extension in latitude of both the EGRET 
[27] and MILAGRO emission [123]. In this case the extension 
of both GeV and multi-TeV emission is determined only by 
the gas distribution in the Galactic disk, which constitutes the 
target for proton-proton interactions. In the leptonic scenario, 
the multi-TeV emission is produced via inverse Compton 
scattering by ~ 100 TeV electrons which can propagate over 
a distance of w 100 pc before being cooled by synchrotron 
and inverse Compton losses [122]. Thus, also in this case 
a narrow latitude distribution of gamma rays is expected, if 
the sources of cosmic ray electrons are concentrated around 
the galactic plane. In this case, however, the extension of the 
gamma-ray emission would be energy dependent, since higher 
energy electrons are cooled faster and can propagate shorter 
distances [122]. Thus, an accurate comparison of the extension 
in latitude of the emission detected at GeV and multi-TeV 
energies is of crucial importance, and future observations with 
improved angular resolution by FERMI and HAWK might 
help in discriminating between the two scenarios. 



The fact that both the Cygnus region and the inner Galaxy 
are characterised by an enhancement in the gas density might 
suggest an hadronic origin for the diffuse multi-TeV emission 
detected by MILAGRO. Moreover, in the case of the Cygnus 
region, the morphology of the multi-TeV emission correlates 
with the CO emission, which traces the gas density [29]. 
Similarly, for the region close to the inner Galaxy, the rather 
narrow latitude profile of the multi TeV-emission, concentrated 
around the dense Galactic disk, might suggest a hadronic 
interpretation of the gamma-ray data, though also a leptonic 
one seems feasible [123]. Since the production of gamma rays 
is accompanied by a corresponding neutrino flux only in the 
hadronic scenario, the search for a diffuse flux of neutrinos 
from the inner Galaxy, which might be detectable by km'^- 
size neutrino detectors, will provide a unique opportunity to 
disentangle the origin of the multi-TeV radiation [124]. 

V. OBSERVATIONS IN THE MULTI-TeV RANGE: 
ARE SNRS COSMIC RAY PEVATRONS? 

The observed CR spectrum extends as a featureless power 
law up to the knee at an energy of a few PeV, where it 
undergoes a significant steepening. This suggests that the 
sources of galactic CRs, whichever they are, must be able to 
accelerate particles up to at least a few PeV. As seen in Sec. 2, 
both flie X-ray observations of SNRs [16], [18], [19] and the 
recent theoretical studies of magnetic field amplifications due 
to CR streaming at shocks [68], [69] seem to suggest the 
presence of a strong magnetic field with intensity ~ 100 fiG 
or more, which is the value needed to reach (or even exceed) 
PeV particle energies in diffusive shock acceleration. 

A decisive and unambiguous indication for the acceleration 
of PeV protons in SNRs can be provided by observations of 
7-rays at energies up to 100 TeV and beyond. Because of the 
Klein-Nishina effect the efficiency of inverse Compton scat- 
tering in this energy band is dramatically reduced. Thus, the 
possible detection of SNRs in gamma rays up to these energies 
could be unambiguously interpreted as an evidence for the 
fact that such objects act as CR PeVatrons [22]. Although the 
potential of the current ground-based instruments for detection 
of such energetic 7-rays is limited, it is expected that the next 
generation arrays of imaging Cherenkov telescopes optimized 
in the multi-TeV energy range will become powerful tools for 
this kind of studies [50], [125]. 

It should be noted that the number of SNRs currently 
bright in > 10 TeV 7-rays is expected to be rather limited. 
Multi-PeV protons can be accelerated only during a relatively 
short period of the SNR evolution, namely, at the end of 
the free-expansion phase/beginning of the Sedov phase, when 
the shock velocity is high enough to allow sufficiently high 
acceleration rate. When the SNR enters the Sedov phase, 
the shock gradually slows down and correspondingly the 
maximum energy of the particles that can be confined within 
the remnant decreases. This determines the escape of the 
most energetic particles from the SNR [26]. Thus, unless our 
theoretical understanding of particle acceleration at SNRs is 
completely wrong, we should expect an energy spectrum of 



CR inside the SNR approaching PeV energies only at the 
beginning of the Sedov phase, typically for a time < 1000 
years. When the SNR enters the Sedov phase, the high energy 
cutoffs in the spectra of both protons and 7-rays gradually 
moves to lower energies, while the highest energy particles 
leave the remnant [26]. This can naturally explain why the 7- 
ray spectrum of the best studied SNR RX J1713.7-3946 above 
10 TeV becomes rather steep with photon index ps 3 [49]. 

The detection of multi-TeV 7-rays generated by the CRs 
that escape the SNR might constitute the smoking gun for 
particle acceleration up to PeV energies [22]. A molecular 
cloud located close to the SNR can provide an effective target 
for production of 7-rays. The highest energy particles few 
PeV) escape the shell first. Moreover, generally they diffuse 
in the interstellar medium faster than low energy particles. 
Therefore they arrive first to the cloud, producing there 7- 
rays (and also neutrinos) with very hard energy spectra. Note 
that an association of SNRs with clouds is naturally expected, 
especially in star forming regions [96] . The duration of 7-ray 
emission in this case is determined by the time of propagation 
of CRs from the SNR to the cloud. Therefore 7-ray emission 
of the cloud lasts much longer than the emission of the SNR 
itself. This makes the detection of delayed 7-ray and neutrino 
signals more probable. The detection of these multi-TeV 7- 
rays from nearby clouds would thus indicate that the nearby 
SNR in the past was acting as an effective PeVatron [22]. 
Both 7-rays and neutrinos are emitted with fluxes detectable 
by currently operating and forthcoming instruments. Since the 
7-ray spectra from clouds are extremely hard, 7-ray telescopes 
operating at very high energies (> 10 TeV) would be the best 
instruments for this kind of study. Remarkably, a detection of 
such emission would not only reveal the acceleration of PeV 
CRs, but also suggest the best targets for neutrino observations, 
since neutrino telescopes are optimized for observations in the 
multi-TeV energy region. 

VI. CONCLUSIONS AND FUTURE PERSPECTIVES 

The recent developments of gamma ray astronomy, espe- 
cially at very high (TeV) energies, constitute an important step 
towards the solution of the problem of the CR origin. The de- 
tection of SNRs in TeV gamma rays by Cherenkov telescopes 
allows us to study several aspects of particle acceleration in 
these objects with unprecedented accuracy (see e.g. [11]). On 
the other hand, the observation of diffuse TeV emission from 
the galactic plane by air shower detectors such as MILAGRO 
[123], when combined with earlier EGRET observations of 
the diffuse GeV emission [27], provide useful information on 
the spatial distribution of CRs in the Galaxy and on their 
spectrum. The combination of pointed observations performed 
with Cherenkov telescopes and the continuous monitoring 
of the whole sky by air shower detectors revealed to be a 
very successfull approach, and the complementarity of these 
two techniques will hopefully be exploited in the future by 
instrument Uke CTA [113], AGIS [114] and HAWC [126]. 

Though both morphological and spectral studies of SNRs 
seem to favor an hadronic origin of the gamma ray emission. 



leptonic models still cannot be considered ruled out. The 
detection of neutrinos from the direction of SNRs would 
unambiguously solve the issue, and prove that SNRs can 
indeed accelerate CR protons [9], [20]-[23]. However, since 
such detections appear to be challenging even for km^-scale 
neutrino telescopes, the search for evidence of CR proton 
acceleration coming form gamma ray observations is manda- 
tory. Besides the issue of the hadronic or leptonic origin of 
the gamma ray emission from SNRs, the problem of the 
(unexpected) high level of isotropy observed from CRs up 
to the energy of the knee stiU persists (e.g. [25]). 

The solution of the long standing problem of the origin of 
galactic CRs might finally come from an improved knowledge 
of the CR spectrum and spatial distribution throughout the 
Galaxy and from a better understanding of several aspects of 
the problem, including acceleration of CRs at shocks, escape 
of CRs from SNRs and their propagation in the Galaxy. 

A. Future perspectives: theory 

The non Unear theory of diffusive shock acceleration gives 
us a satisfactory description of how particles are accelerated 
at shocks and has been successfully applied to model the 
multiwavelength emission from SNRs (see e.g. [75], [89]). 
Despite the success of the theory, some issues still need further 
investigation. In particular, new aspects of the basic physics 
of the CR-driven magnetic field ampUfication at shocks have 
been understood [69], but a fully self consistent treatment of 
the problem is still missing. A better understanding of this 
mechanism might provide a solid theoretical interpretation for 
the very high values of the magnetic field inferred from X-ray 
observations of SNRs [16]-[19] and will allow us to estimate 
the particle acceleration rate and thus the maximum particle 
energy achievable within this theoretical framework [69]. 

Another missing piece of information is the way in which 
accelerated CRs leave the SNR. While protons are believed 
to be released gradually during the whole Sedov phase as the 
SNR shock slows down [26], electrons are trapped within the 
remnant due to severe synchrotron losses and are probably 
released in the interstellar medium solely during the late 
phase of the SNR evolution. Surprisingly, this aspect of the 
problem has not yet received as much theoretical attention 
as it deserves. Studies in this direction will allow us to 
estimate the spectrum of CR protons and electrons injected 
in the interstellar medium and will possibly tell us something 
about how the proton to electron ratio varies during the SNR 
evolution. 

After being released into the interstellar medium, CRs 
propagate in the galactic magnetic field. The propagation of 
CRs is generally investigated by means of numerical codes 
that solve the transport equation given the diffusion coefficient 
and possibly the structure of the galactic magnetic field. The 
most popular amongst these codes is GALPROP [36], which 
provides as an output also the diffuse garmna ray emission 
from CR interactions in the interstellar gas and radiation field 
[120]. The main limitation of such approaches is that they 
normally assume a continuous and spatially smooth injection 



of particles along the galactic plane and, being normalized 
to local CR data, cannot account for localized excesses of 
CR sources. These limitations became evident when attempts 
to use GALPROP to model the diffuse emission detected by 
MILAGRO from the Cygnus region failed to fit observations 
by about an order of magnitude [29]. Thus, though GALPROP- 
like approaches can predict global properties of the diffuse 
gamma ray emission from the Galaxy, they cannot be applied 
to model local excesses of the diffuse gamma ray emission. 

One possible way to improve the studies of CR propagation 
and of the related diffuse emission is to include in numerical 
codes the stochasticity of SNR explosions, which would lead 
to the appearance of localized (both in time and space) 
excesses in the CR density close to their sources [22], [106], 
[127]. The expected association between supernovae and dense 
gas regions that provide a thick target for CR hadronic 
interactions might further enhance the resulting gamma ray 
emission. Obviously, this approach, being statistical in nature, 
would only allow to infer average properties of the gamma ray 
emission from the Galaxy, such as, for example, the number of 
molecular clouds irradiated by CRs detectable in gamma rays 
as unidentified sources [107], or the amplitude and the spatial 
scale of the fluctuations in the diffuse gamma ray emission. 
These predictions would provide clear observational tests for 
propagation models and for the SNR hypothesis for CR origin. 

Finally, self consistent propagation models where CRs 
propagate away from their sources due to scattering on the 
magnetic turbulence they themselves produce [128] need to 
be further developed and, together with the knowledge of the 
CR spectrum injected by SNRs in the interstellar medium, 
might clarify the issue of the isotropy observed in the arrival 
direction of CRs up to PeV energies. 

B. Future perspectives: observations 

Future observations in the gamma ray domain are expected 
to drive the progresses in the field. Comprehensive reviews 
of the expected performances and scientific outputs from next 
generation gamma ray instruments can be found in Refs. [11], 
[125], [129]. Here the expected contribution from these obser- 
vations to the problem of CR origin is summarized. 

First of all, it has to be noted that observations in different 
energy regions will provide insights into different aspect of 
the problem. Following previous classifications, we divide the 
gamma ray spectrum into three regions: i) the GeV domain, 
which spans the energy spectrum in the interval 100 MeV - 
100 GeV, currendy covered by the FERMI sateUite; ii) the TeV 
domain, where Cherenkov telescopes and air shower detectors 
operate, w 100 GeV - 10 TeV, and Hi) the multi-TeV domain, 
approximatively > 10 TeV, which will be explored by next 
generation of instruments. 

1 ) GeV Energy Domain: FERMI: FERMI is now scanning 
the whole sky at GeV energies with unprecedented sensitivity. 
The recent detection of the SNR IC443 by AGILE suggests 
that FERMI, with its better sensitivity will possibly detect 
several of these objects. A detection of SNRs at GeV energies 
is of crucial importance to discriminate amongst hadronic and 



leptonic origin of the gamma ray emission. Moreover, FERMI 
will measure the diffuse emission from the galactic plane, 
and thanks to its improved angular resolution is expected 
to reveal structures in the diffuse emission which reflect the 
underlying dumpiness of the interstellar medium. The possible 
detection of individual molecular clouds will allow us to probe 
the CR spectrum in different parts of the Galaxy, providing 
stringent constraints to CR propagation models and revealing 
the presence of localized excesses of CR sources. 

2) TeV Energy Domain: how many sources?: Next gen- 
eration of instruments operating in the TeV domain are ex- 
pected to increase significantly the number of sources detected 
along the galactic plane. The improved angular resolution of 
Cherenkov telescopes will allow better morphological studies 
of extended galactic sources, useful to search for correlations 
between gamma ray emission and gas density and/or X- 
ray emission. The presence or absence of such correlations 
will help in discriminating between hadronic and leptonic 
models for the gamma ray emission. Air shower detectors will 
complement such studies by mapping the extended emission 
from the whole Galaxy. 

3) Multi-TeV Energy Domain: searching for Cosmic Ray 
PeVatrons: Most of the galactic sources detected so far by 
Cherenkov telescopes do not show evidence for a high energy 
cutoff in their spectra [130]. For this reason, an extension 
of the observed energy range in the multi-TeV domain is 
desirable, since it will allow us to search for cutoffs in such 
spectra, thus probing the most extreme particle accelerators in 
the Galaxy. Moreover, the sources of CRs with energies up 
to the knee are expected to exhibit a cutoff in this domain 
of the gamma ray spectrum. This implies that observations in 
the multi TeV region might finally reveal the nature of the CR 
PeVatrons. 
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